Erratum: Nodal Cooper-pair stabilized phase dynamics in granular (i-wave 

superconductors 



Yogcsh N. Joglekar, A. H. Castro Neto, and Alexander V. Balatsky 



Recently, we discovered a mistake in our Letter 111 on phase dynamics in granular (i-wave superconductors, which 
we corrected in an Erratum published last week |2j ■ This is an expanded version of the Erratum, containing details 
and supporting numerical results. 

In our Letter, we accounted for weak disorder in d-wave superconducting state by adding lifetimes to the nodal 
quasiparticles E^ — > E^ — i^y (7 > 0) when evaluating thermal Green's functions G(k, r) and T(k, r). This naive 
substitution in thermal Green's functions leads to Gdis ~ Gexp(i7r) and i*dis ~ Texp(z7T). However, the exact 
expressions for thermal Green's functions at temperature T are given by 
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where = (e^ — yu) 2 + A 2 , is the quasiparticle dispersion in the superconducting state, Ak = Ao cos 20k is the 
(i-wave superconducting gap, u>„ = Tn(n + 1/2) are fermionic Matsubara frequencies, and we have characterized the 
self-energy due to disorder as $sY,(iu) n ) — — 7sign(u; n ). 

The thermal Green's functions, Eqs.Q) and are purely real and satisfy the following relations 

G(k,r) = -G(k,-r) = +G(k, 1/T - r), (3) 
F(k,r) = +F(k,-r) = -F(k,l/T-r), (4) 

where -1/T < r < 1/T. Consequently, they also satisfy G(l/2T-r) = G(l/2T+r) and F(1/2T-t) = -F(l/2T+r) 
where |r| < 1/2T. The Matsubara frequency sums above can be re-expressed in terms of sum-over-residues, and an 
integral over a branch cut which cannot be performed analytically. In the T = limit, the sums reduce to 
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These integrals can be expressed in terms of Sine (I s ) and Cosine (I c ) integrals 

G(k,r) = -i-[7 s ( r , 7 -^ k )+/ s (r, 7 + ^ k )], (7) 

F(k,r) = --^|-[/ c (T,7-*£k)-/ c (T,7 + ^k)]. (8) 

but the results, Eqs.Q and JSJ), are not particularly illuminating. They can be expressed in terms of elementary 
functions only in the clean limit 7 — ► 0, in which case they reduce to standard analytical results. 

It is straightforward to evaluate the frequency sums in Eqs.Q and J3J numerically. Figure shows G(r) and F(t) 
for — 1/T < t < 1/T for different disorder strengths. The Reader can see that the numerical results obey all symmetry 
properties listed in the last paragraph. At low temperatures, the long-time behavior of the disordered thermal Green's 
functions is shown in Fig|21 It is clear that, for T <C Ao, 7 and tT <C 1, both G(r) and F(t) decay with disorder. We 
can approximate them as Gdi S {r) = G(r) cxp(— cig7|t|) and Fdi s (r) — F(t) cxp(— of 7 |t|), where ac, o_f ~ 0(1) (See 
Fig. © . 

Therefore, although Eqs.(9) and (12) in [l( are correct, our results ad%s{r) = cos(27t)«(t), fidisir) = j3(r) are 
incorrect. Instead, we find that the disorder suppresses the power-law correlations in Eqs.(9) and (12) inQ, 

a dis (T) ~ {G L G R + G* R G* L ) = e- 2a ^a(T), (9) 
Mr) ~ [F L F* R + F R F* L ) = e- 2a ^/3(r). (10) 
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Since the disorder suppresses the long-range phase correlations, it is necessary to exercise caution in using gradient 
expansion to map the nonlocal imaginary-time phase action (Eq.(4) in on to a local Lagrangian for a Joseph- 
son junction (Eq.(17) in Q). In particular, the (apparently) divergent contribution to capacitance, which arises 
from momentum-conserving tunneling term (Eq.(16) in is incorrect. Instead the junction capacitance is deter- 
mined by the point-contact tunneling term, C ~ A)Ao ~ |To| 2 /Ao, and the Josephson energy is determined by the 
momentum-conserving tunneling term. Subsequently, our prediction that the Josephson plasmon frequency is strongly 
renormalized due to long-range phase correlations and is not sensitive to disorder, ui p ~ y/AoT, Eq.(20) in Q is not 
valid. The correct dispersion is uj p ~ Ao and it is weakly renormalized due to disorder. 

Thus, our results on robustness of d-wave Josephson junctions against disorder and the strong renormalization of the 
Josephson plasmon frequency are invalid. The clean-limit results in our paper, such as power-law phase-correlations 
and the subsequent sub-ohmic dissipation, continue to remain valid. 

We wish to thank Professor Dmitri Khveshchenko for his questions which led to the discovery of our mistake. 
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FIG. 1: Normal and anamolous propagators, G(r) and F(t), calculated numerically for different disorder strengths 7 by using 
Eqs. Q and @. The disordered Green's functions satisfy the symmetry properties in Eqs.JSJ and Q. 




FIG. 2: G and F functions for T <C E, 7. It is clear that both decay with increasing disorder. The kink in G(r) at small r is 
due to the discontinuity in the single-fermion propagator at r = 0. 
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FIG. 3: Dependence of G(r) and F(t) on disorder, characterized by using the fit Gnu = Gexp(— cig7|t|) and F&i a = 
_Fexp(— a_F7|r|), for E = 10. In the present case, ac = 0.6 and ap = 0.5. The intercept on the j/-axis is — E = —10 
and is consistent with the analytical result for low temperatures. 
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